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ABSTRACT 

Ph Context. As is the case for all techniques involved in the research for exoplanets, direct imaging has to take into account the probability 

PJ of so-called astrophysical /a/^e positives, which are phenomena that mimic the signature of objects we are seeking. 

Aims. In this work we aim to present a case of a false positive found during a direct imaging survey conducted with VLT/NACO. A 
promising exoplanet candidate was detected around the K2-type star HD 8049 in July 2010. Its contrast of AH = 7.05 at 1.57 arcsec 
allowed us to guess the presence of a 35 Mjup companion at 50 projected AU, for the nominal system age and heliocentric distance. 
Q Methods. To check whether it was gravitationally bound to the host star, as opposed to an unrelated background object, we re-observed 

^ the system one year later and concluded a high probability of a bound system. We also used radial velocity measurements of the host 

"^ star, spanning a time range of ~ 30 yr, to constrain the companion's mass and orbital properties, as well as to probe the host star's 

jrt spectral age indicators and general spectral energy distribution. We also obtained for the companion (/-band imaging with EFOSC 

I I and NIR spectroscopy. 

Results. Combining all these information we conclude that the companion of HD 8049 is a white dwarf (WD), with temperature 
^~^ T^ff = 18800 ± 2100 K and mass Mwd = 0.56 ± 0.08 Mq. The significant radial velocity trend coupled with the imaging data indicate 

^ that the most probable orbit has a semimajor axis of about 50 AU. The discrepancy between the age indicators suggested against a 

^^ bona-fide young star. The moderately large level of chromospheric activity and fast rotation, mimicking the properties of a young 

C^ star, might be induced by the exchange of mass with the progenitor of the WD. This example demonstrates some of the challenges in 

^~H determining accurate age estimates and identifications of faint companions. 

• Key words. Techniques: high angular resolution - spectroscopic - Stars: white dwarfs - brown dwarfs - individual: HD 8049 - binaries: 
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. . 1 . Introduction e.g. Queloz et al. 2001; Lovis et al. 201 1). A number of contro- 

^ versial cases have been reported in the literature (e.g. Setiawan 

- In the search for extrasolar planets, we are frequently faced ^^ ^j 20O8; Huelamo et al. 2008; Hernan-Obispo et al. 2010; 

with objects or physical phenomena that mimick the signature pigueira et al. 2010). Simultaneous monitoring of activity indi- 

^ of a substellar object — the so-called astrophysical false pos- ^^^^^^^ ^^ec\ia\ line profile changes and photometric variations 

itives. The most commonly used techniques in the search for ^^^ ^jj^^ one to recognize the true origin of the radial veloc- 

exoplanets are the radial velocity, transit, TTV (Transit Timing ^^y variations and to at least partially account for them (see e.g. 

Variation), microlensing and direct imaging techmques. A large goiggg ^^ ^1 201 1 ; Lanza et al. 2010; Dumusque et al. 2012). 
fraction of the planets discovered so far has been detected 

through the radial velocity technique. Rotational modulations of Transit searches are even more heavily plagued by false pos- 
star spots and magnetic activity cycles might produce spurious, itives. Various kinds of eclipsing binaries can mimick a transit- 
periodic, radial velocity signals that mimic planet signatures (see ing planet photometric curve. Such false positives can often out- 
number the true exoplanet detections (Brown 2003). The phe- 






X 



* Based on observations collected at La SiUa and Paranal nomenon of blended eclipsing binaries is the most severe, as 

Observatory, ESO (Chile): Programs 184.C-0567 (NAOS-CONICA), described in O'Donovan et al. (2006). To help mitigate these 

184.D-1151 (EFOSC), 60.A-9036 (HARPS), 089.C-0665 (SINFONI). adverse effects, various observational diagnostics, based on the 

Based on the All Sky Automated Survey (AS AS) photometric data. photometric light curves or on additional spectroscopic or pho- 



Zurlo et al.: Astrophysical false positive in the NACO 



tometric observations, have been developed. Statistical tools can 
also be exploited to explore all possible scenarios of binary sys- 
tems and the probability of their occurrence (e.g. BLENDER 
code, Torres et al. 2011, PASTIS code, Diaz et al., in prepara- 
tion). 

Timing variations in eclipsing binaries might be signatures 
of circumbinary planets, but alternative explanations cannot be 
firmly ruled out (Zorotovic & Schreiber 2012). A recent case of 
TTVs induced by stellar activity (star spot features modulating 
the brightness of the host) instead of the gravitational influence 
of a planet is described in Barros et al. (2013). Finally, a case of 
a possible microlensing planet detection was recently revealed 
to be a false positive, resulting from the variability of the lensed 
star (Gould et al. 2012). 

Concerning the direct imaging technique, the subject of this 
work, it may occur that the detected candidate companion is not 
gravitationally bound to the host star Detection of unbound ob- 
jects is frequent since direct imaging surveys typically have bet- 
ter sensitivities to to planets far from the host star, where a higher 
contrast can be reached (see e.g. Chauvin et al. 2010; Vigan et al. 
2012). If the follow-up observation is close in time to the first 
epoch and the displacement of the star on the sky due to proper 
motion is small, it may be challenging to determine with a high 
probability whether an object is bound. One example is given by 
Metchev & Hillenbrand (2009) where an object that was origi- 
nally thought to be compatible with being bound, based on the 
proper motion test, was determined through spectroscopy to be 
a background M-star Moreover, the case of the false positive 
object around IM Lupi (Mawet et al. 2012) teaches that, when 
the distance determination is complex, it is possible to miscalcu- 
late the relative certainties of a bound versus unbound interpre- 
tation. Secondly, even after a candidate companion is shown to 
be bound, the system might prove to be a stellar binary, which is 
the case of the object of our study. 

Finally, in order to meaningfully characterize the compan- 
ion, the system age must be accurately determined, as the lumi- 
nosity of substellar objects depends strongly on the age of the 
system (younger age corresponds to a brighter companion mag- 
nitude, see e.g. Baraffe et al. 2003; Madhusudhan et al. 2011); 
spurious young stars might therefore lead to poorly classified 
companions. In this paper we present the case of such a spuri- 
ous young star, along with an originally misidentified substel- 
lar candidate, found in the NACO Large Program (NACO-LP'; 
Chauvin 2010). 

The object of our study, HD 8049, was observed in July 
2010. The target is a solar neighborhood K2-star, with a distance 
of 33.6 pc. The companion found from the direct imaging obser- 
vations had an //-band magnitude consistent with a substellar 
mass, based on the nominal age of the central star. However, af- 
ter a deeper analysis, presented in this paper, it is concluded to 
be a white dwarf (WD). The case is similar to that of the ob- 
ject orbiting Gliese 86 (see Mugrauer & Neuhauser 2005), ini- 
tially considered a brown dwarf by Els et al. (2001). Other cases 
of close WD companions around HD 147513 and HD 27442 are 
described in Porto de Mello & da Silva (1997) and Chauvin et al. 
(2007) respectively. 

Sometimes, as in the case presented here, a false positive 
turns out to be an interesting object: among known nearby WD 
this target is in fact one of the hottest. Because it is part of a close 
binary system, it yields insights into how those types of object 
interact with their companions. 



Table 1. Resuming table of the observations of HD 8049. 



Instrument 


Date 


Type of Data 


CORAVEL 


Iunel4UT, 1983 


1 RV p.,V spectrum 




December 31 UT, 1983 


1 RV p.,V spectrum 




November 30 UT, 1985 


1 RV p.,V spectrum 


CORALIE(98) 


from December 26 UT, 2000 
to January 21 UT, 2007 


43 RV points 


HARPS 


September 15 UT, 2004 


1 RV p.,V spectrum 


CORALIE(07) 


from July 3 UT 2007 
to July 8 UT, 2012 


20 RV points 


NACO 


July 30 UT 2010 


//-band imaging 




July 28 UT 2011 


//-band imaging 


EFOSC 


December 21 UT 2011 


(/-band imaging, 
V spectrum 


SINFONI 


August 28 UT, 2012 


JHK spectra 



' ESO program 184.C-0567, PI J.-L. Beuzit, "Probing the 
Occurrence of Exoplanets and Brown Dwarfs at Wide Orbits". 



This work is structured as follows: we present the observa- 
tions in Section 2; we show the companion characteristics in 
Section 3; we summarize the stellar properties and how they are 
influenced by the companion in Section 4; the properties of the 
binary system are discussed in Section 5; and conclusions are 
drawn in Section 6. 



2. Observations and data reduction 

HD 8049 was observed with several instruments and techniques. 
All data we collected for this analysis are listed in Table 1 along 
with the related instrumentation, date of observation and tech- 
nique. 

2.1. NACO observations 

HD 8049 was observed as part of the NACO Large Program 
(NACO-LP) imaging survey. We used VLT/NACO (Lenzen et al. 
2003; Rousset et al. 2003) high-contrast Adaptive Optics (AO) 
system with the //-band (1.33-1.99 fim) and the S13 camera, 
which provides a spatial sampling of ~13 mas/pixel and a field 
of view (FoV) of 14" x 14". The target was observed in pupil- 
tracking mode to allow implementation of angular differential 
imaging (ADI, Marois et al. 2006) data analysis techniques. The 
observing sequence consisted of a series of unsaturated images 
of the point spread function (PSF), which served as reference 
for differential photometry and astrometry, and a series of deep 
saturated exposures optimized for the detection of faint com- 
panions. Data reduction for saturated and unsaturated sequences 
followed standard procedures (flat-fielding, sky subtraction, bad 
pixel coiTection). The saturated sequence was then analysed with 
two methods: (1) the LOCI algorithm (Lafreniere et al. 2007) 
with a separation criterion A^^ - 0.75 FWHM, and (2) a spatial 
filtering of all images in a 5x5 FWHM box followed by de- 
rotation and median-combination. Differential astrometry was 
performed using fitting of a 2D Gaussian profile on the primary 
star and on the point sources detected in the reduced images. 

The target was first observed on July 30 UT, 2010. A 
faint point source (hereafter HD 8049 B) at a separation of p = 
1.566 ± 0.006" and with a position angle of 0= 118.4 + 0.2 deg 
was detected with a magnitude difference AH = 7.1 ± 0.1 mag 
(see Figure 1). Follow-up observations on July 28 UT, 2011 re- 
vealed that the point source was co-moving with the primary. 
Thanks to the high proper motion of the central star (//„ = 
65.99 ±1.18 mas/yr and fig = -240.99 + 0.98 mas/yr), a x^ 
probability test on Aa and AS with respect to the star at two 
epochs rejected the possibility (at 99% certainty) that the object 
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was a background source. The measured contrast, along with our 
original age estimate and the measured astrometry, led us to con- 
clude the presence of a 35 Mjup companion at 50 AU projected "^^^ 
separation. 




Fig. 1. VLT/NACO, i/-band image (1.65 yum) of the companion. 
The frame was generated from the median combination of pro- 
cessed images having undergone spatial filtering (5x5 A/ D box) 
and de-rotation. The central core of the PSF is saturated and has 
been masked in this image. 




Fig. 2. Radial velocity of the star measured with four instru- 
ments: CORAVEL (green), CORALIE98 (red), CORALIE07 
(blue) and HARPS (black). The solid black line is a least-squares 
fit of a quadratic RV trend. While the drift information alone is 
consistent with either a BD or WD interpretation, the inclusion 
of the separation information, as gathered from NACO, allows 
us to rule out the BD possibility. 



2.2. CORAVEL, CORALIE and HARPS high resolution 
spectroscopy 

HD 8049 was observed as part of a long term extra-solar planet 
search program started in June 1998 with the CORALIE high- 
resolution fibre-fed echelle spectrograph mounted on the 1.2 m 
Euler Swiss telescope at La Silla (ESO, Chile). In addition to the 
CORALIE data, we obtained a spectrum with HARPS. Finally, 
we used the CORAVEL database (Baranne et al. 1996) and ex- 
tracted 3 radial velocity (RV) observations with a 300 m/s preci- 
sion (Baranne et al. 1996). This enabled us to have a total time- 
span of nearly 29 years. 

For these three instruments, radial velocities were computed 
by cross-correlating the observed stellar spectra with a binary 
mask, whose non-zero zones correspond to the theoretical posi- 
tions and widths of stellar absorption lines at zero velocity. For 
HARPS and CORALIE, the instrumental velocity drifts were 
monitored and corrected using the "simultaneous thorium ref- 
erencing technique" with dual fibres (more details in Baranne 
et al. 1996), whereas for CORAVEL a standard radial velocity 
star (HD 168454) was used for calibration. Radial velocity points 
are shown in Figure 2. The trend in the radial velocity points is 
likely produced by the influence of a more massive object than a 
brown dwarf. Additionally, activity indicators (S index, log ^hk) 
and the projected rotational velocity were also measured on the 
CORALIE and HARPS spectra. 

2.3. NTT/EFOSC 

We obtained 10 images of HD 8049 on December 21 UT, 201 1, 
with EF0SC2 at the NTT (ESO, Chile). Each frame was taken 
with the Johnson U filter, a 0.2 s exposure to avoid saturation, 
2x2 binning and seeing of 0.8". We reduced the {/-band im- 
ages with standard IRAF- tasks. We fitted the two objects with 



Moffat profiles and performed aperture photometry on each star 
after subtracting contributed flux from the neighbor star The 
resolved binary system had a magnitude difference of AU = 
1.714 ± 0.023 mag, as shown in Figure 3. The difference was 
calculated taking into account the correction for the U - B color 
index of the two stars using equation (9) of the EFOSC2 man- 
ual \ We used for the K star the value listed in Table 15.7 of 
Allen (1977) and for the WD we considered the colors in the 
range of effective temperature, Tf-g, from 15000 K to 25000 K 
(see Sec 3.2). The measured values of the projected separation, 
p= 1.56 + 0.02", and the position angle, = 118±0.5deg, were 
compatible with NACO measurements. 




Fig. 3. [/-band EFOSC image of the binary system. HD 8049 B, 
in the south east side, is 1 .7 mag fainter. 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 



for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
^ EF0SC2 manual, version 3.6. ESO user facilities. 
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During the same night, we performed spectroscopy of the 
system. The n.ll grism and the 1" slit were used, providing a 
resolution R ~ 360 over the range between 3500-7400 A. The 
slit was oriented along the two stars (i.e. PA = 121 deg). We 
obtained two frames of 8 s exposure. 

The spectral reduction was done using standard IRAF tasks. 
The raw data were bias subtracted and flat-field corrected. The 
spectra were extracted with the IRAF task apall, wavelength cal- 
ibrated using arc lamp exposures, and cross-correlated with the 
sky lines. The relative flux calibration was performed by deriv- 
ing the sensitivity curve for the instrument set-up from the spec- 
troscopic standard star L745-46A (Hamuy et al. 1994). 

We found that a low S/N signature of HD 8049 B was present 
in the blue part of the spectrum; its contribution is completely 
blended with that of the primary star. 

2.4. VLT/SINFONI 

The companion was further characterized on August 28 UT, 
2012 using the Spectrograph for INtegral Field Observations in 
the Near Infrared (SINFONI, see Eisenhauer et al. 2003; Bonnet 
et al. 2004b,a) mounted at the VLT/UT4. SINFONI provides 
AO assisted integral field spectroscopy in the near-infrared. It is 
composed by the Spectrograph for Infrared Faint Field Imaging 
(SPIFFI) fed by a modified version of the Multi-Applications 
Curvature Adaptive Optics system (MACAO, see Bonnet et al. 
2003). We used the gratings J and H + K to cover the 1.1- 
1.4 and 1.45-2.45 |/m ranges at a spectral resolution of ~2000 
and ~1500 respectively. The smallest plate scale (12.5 mas) was 
used, enabling Nyquist sampling of the PSF The primary served 
as a reference for the AO wavefront sensor and was kept outside 
the instrument field of view (0.8" x 0.8"). Each sequence was 
composed of 12 acquisitions on the source with small dithering 
and one acquisition on the sky. Immediately after observing our 
science targets, we observed the standard star HIP 17280 (B5V) 
following an AB pattern to correct our spectra for telluric fea- 
tures. HD 8049 and the standard were both observed under poor 
atmospheric conditions with a seeing of 1.3-1.5", although the 
AO loop remained stable, and airmass was lower than 1.1. 

We homogeneously reduced our dataset with the ESO data 
reduction pipeline version 1.9.8 (Abuter et al. 2006). The 
pipeline successfully carried out bad pixels detection, flat field- 
ing, distortion coefficient computation, and wavelength calibra- 
tion using calibration frames acquired the day after the ob- 
servations. Individual datacubes were reconstructed from sky- 
subtracted object frames and merged into a final mosaicked dat- 
acube. Datacubes of telluric standard stars were obtained in a 
similar way. Finally, we integrated the source flux in each final 
datacube over selected apertures to optimize the signal to noise 
ratio without introducing differential flux losses. The spectra of 
the standard star were corrected from their intrinsic features and 
divided by a black body at the appropriate temperature (Pickles 
1998). The SINFONI near-infrared spectrum of the object is 
shown in Figure 4. 

3. Nature and properties of HD 8049 B 

3.1. Early discovery 

HD 8049 was identified as a binary system by van den Bos 
(1929) in 1928. A projected separationp = 3.19" and a position 
angle = 153.1 deg were measured at epoch 1928.81. We cor- 
rected for the precession effect as listed in Heintz (2000) and find 
6 - 153.6 deg. The original paper did not provide any error bar 



on the measurement, just a note that there were good observing 
conditions. We assumed an error of 100 mas on Aa and A6. The 
magnitude resulting from the visual observations was 13.5 mag 
(rounded to 0.5 mag) and a bluish color was reported; the as- 
sumed error is 0.5 mag. This is the only available measurement 
of the pair listed in the Washington Double Star (WDS) Catalog 
(Mason et al. 2001) before our observations. This position is pos- 
sibly compatible with the one measured in our EFOSC [/-band 
images and with the NACO results, once the time lapse is con- 
sidered. Table 2 summarizes the available astrometry, including 
the two epochs from NACO. 

From NACO observations spanning over 1 yr, we derive 
a relative motion of 2.7 mas/yr in a and 30.4 mas/yr in 6. 
Doing the same for 1928 vs 2010 observations we get a mo- 
tion of 0.4 mas/yr in a and 26 mas/yr in 6. These determinations 
are obviously much smaller than the proper motion of the star 
ifia = 65 .99 ±1.1 8 mas/yr and ;U5 = -240.99 ±0.98 mas/yr), sup- 
porting the physical association of the faint companion. The mo- 
tions with respect to the central star derived from 1928 to 2010 
and from 2010 to 2011 are consistent within their errors, mak- 
ing viable the hypothesis that the companion detected in 1928 in 
the optical range is the same one we found with NACO, and that 
we are indeed tracking the orbital motion of HD 8049 B around 
HD 8049 A. The characterization of the orbital properties of the 
binary is discussed in Section 5.1. 



Table 2. Astrometry of candidate(s) around HD 8049 



Epoch 


Aa 


AS 


Reference 




(mas) 


(mas) 




1928.81 


1417 ± 100 


-2858 + 100 


van den Bos (1929) 


2010.58 


1378 ± 7 


-744 ± 7 


NACO 


2011.58 


1382 ± 10 


-705 ± 10 


NACO 


2011.97 


1373 ± 24 


-734 ± 24 


EFOSC 



3.2. HD 8049 B, a new r^earby WD 

Considering the combined optical and NIR photometry of 
HD 8049 B {U, V and H bands) from EFOSC, the WDS Catalog 
and NACO observations, the hypothesis of a sub-stellar com- 
panion is unambiguously excluded. The combined visible and 
NIR photometry indicate HD 8049 B is most likely a WD. The 
SINFONI spectrum indeed confirms this assumption as the ob- 
served spectral energy distribution shows a hot and dense atmo- 
sphere that can be associated with that of a WD. The NIR spec- 
trum is rather featureless as we do not detect any hydrogen or 
helium lines over the whole 1.0-2.5 ;um spectral range. Despite 
the low-SNR, it could indicate a featureless Helium-rich nature 
(DB WD) or a Hydrogen-rich white dwarf (DA) with a rela- 
tively hot atmosphere and therefore less contrasted Hydrogen 
lines (see Dobbie et al. 2005). Visible band spectroscopy should 
enable a better characterization of the WD nature of HD 8049 B, 
but unfortunately the low quality EFOSC spectrum does not pro- 
vide such additional information due to the blend with the much 
brighter primary spectrum. 

To further constrain the physical properties of HD 8049 B, 
we used the emitted fluxes in the FUV and NUV, listed in the 
GALEX Catalog (Martin et al. 2005). In the ultraviolet, we ap- 
plied the correction for the chromospheric emission of the pri- 
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Fig. 4. SINFONI spectrum of HD 8049 B in 7-band (left) and HK-hand (right). 



mary as desciibed in Sect. 5.2. We refer to that section for an 
exhaustive summary of the photometry for the two objects. 

As a preliminary analysis of the object properties, we cre- 
ated empirical sequences using the catalog of nearby WDs by 
Giammichele et al. (2012). When required photometry was un- 
available there, supplemented it with available GALEX FUV 
and NUV magnitudes, and 2MASS J, H and Ks magnitudes. 
The final sample consisted of 107 nearby (< 51 pc) WDs: 22 
with FUV magnitudes, 18 with NUV magnitudes and 84 with 
J, H, Ks magnitudes calibrated by Giammichele et al. (2012). 

The catalog of observed WDs listed mass, Tgff, \ogg (loga- 
rithm of the surface gravity) and fcooh all measured from spec- 
troscopy. As a check, we overplotted these object parameters 
with the expected values from theoretical sequences of Vennes 
et al. (2011), and confirmed that they perfectly fit, as shown in 
Figure 5. 

We noticed that the FUV - NUV color of HD 8049 B did 
not fit with theoretical nor empirical models. The NUV magni- 
tude is more affected by the chromosphere of the K star, but if 
we considered the coiTected value, the discrepancies were even 
greater. 

In order to fit our photometry measurements with the SED 
(Spectral Energy Distribution) of a white dwarf object, we used 
the MCMC (Markov Chain Monte Cario) based code PASTIS 
(Diaz et al., in preparation). To extract characteristics of the 
WD atmosphere, we used models by Koester (2010); these mod- 
els predict the mean intensity of the stellar disk and depend 
on \ogg, Teff and the atmosphere composition; the latter is as- 
sumed to be of pure H in this case. To calculate the observed 
flux measured from our viewpoint on Earth, we multiplied the 
intensity by the solid angle of the star, calculated using the dis- 
tance and the radius of the WD. The latter was obtained from the 
reff-logg-fcooi-mass relation from Holberg & Bergeron (2006); 
Kowalski & Saumon (2006); Tremblay et al. (2011); Bergeron 
et al. (201 1)'* (used later on to also compute the derived param- 
eters). 

The synthetic spectrum was then corrected for interstellar ex- 
tinction using the relations of Fitzpatrick (1999) with Ry - 3.1, 
and an assumed reddening value of ^(b-v) = 0.0106, as listed in 
the GALEX catalog. Finally, the spectrum was integrated in the 
bandpass listed in Table 4 as described in Bayo et al. (2008). We 
did not consider the NUV point for the reason listed previously. 



We also omitted the U point because it can be affected by the 
cut-off of the atmosphere; the V point was rejected because the 
observation was made by eye and there is not an associated eiTor. 

The SEDs of the system and its components are shown in 
Figure 6. The fitted parameters are the distance, \ogg and T^-ff. 
As priors, we assumed a normal distribution with a mean dis- 
tance of 33.6 pc and a standard deviation of 1.4 pc (based on the 
Hipparcos parallax of 29.79 ±1.21 mas, van Leeuwen (2007); 
for logg we used a uniform distribution in the range 7.0-9.0; for 
Teff we used a uniform distribution in the range 10000-40000 K. 

We ran 50 chains of 5 x 10^ steps starting from random 
values of the priors. After removal of the burning phase (steps 
before convergence), we merged all the chains, taking into ac- 
count the correlation lengths. The results of the fitted parameters 
are logg = 7.86 ± 0.16 and T^g = 18800 + 2100 K. The dis- 
tance remained equal to the prior. The derived parameters were 
MwD = 0.56 + 0.08 Mo, Rwd = 0.0145 ±0.0016 Rq, and cooling 
time fcooi - 63^j7 Myr (see Figure 6). Note that the errors of the 
values are the 68.3% confidence interval of the fit, and they do 
not take into account possible model errors. 

The cooling time calculated using the empirical sequences is 
98;^^g Myr (see Fig. 7), in agreement with the SED fit. To esti- 
mate the age of the system we must add the main sequence life- 
time to the cooling age of the WD. This is discussed in Sec. 5.3. 



4. Nature and properties of HD8049 A 

We exploited the spectroscopic data described in Sec. 2.2 and ad- 
ditional data from public archives with the aim of improving the 
age determination and other stellar parameters of the primary. 



4.1. Rotation period 

HD 8049 was observed in the All Sky Automated Survey 
(ASAS^* Variable star catalog; Pojmanski 2002). Archive obser- 
vations span from November 20 UT, 2000 to October 28 UT, 
2009. After removing outliers and averaging observations col- 
lected within short time intervals (less than 2 hr), we were left 
with 583 data points for the rotation period search. The average 
photometric accuracy was 0.02 mag. The brightest V magnitude 
was V = 8.68 ± 0.02 mag. To minimize the effects of active re- 



http://www.astro.umontreal.ca/~bergeron/CoolingModels 



http://www.astrouw.edu.pl/asas/ 
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Fig. 5. Absolute magnitude in different bands (from top-left FUV, NUV, V and H) versus effective temperature for white dwarf 
models of Vennes et al. (2011) (black lines) and sample a of nearby dwarfs collected by Giammichele et al. (2012). Stars represent 
objects with all magnitudes available and diamonds represent the others. Colors indicate the mass from the lowest (blue) to the 
highest (red ones). The continuous blue horizontal line represents the magnitude of HD 8049 B. The corresponding error bars are 
plotted as dashed blue lines. The plots show that the objects of the Giammichele et al. (2012) catalog are not peculiar and are well 
described by the theoretical models of Vennes et al. (201 1). 



gion growth and decay on the periodograms, we sectioned the 
data time series into 14 light curves. The light curve amplitude 
was variable, with values up to AV = 0.04 mag. The Lomb- 
Scargle analysis (Scargle 1982) was carried out on the complete 
time series data as well as on each of the 14 light curves. We 
found a rotation period of P = 8.3 ±0.1 days in the complete 
time series data, and the same period in 4 out of 14 time inter- 
vals (with FAP, False Alarm Probability, < 1%). The light curve 
amplitude was larger than the photometric accuracy cr in 13 out 
of 14 time segments. 

HD 8049 was also observed by SuperWASP (Pollacco et al. 
2006; Butters et al. 2010). The public archive permitted us to 
retrieve data for the 2006 and 2007 seasons^\ We carried out an 
analysis similar to that described for the ASAS data. After out- 
lier removal, averaging, and data series sectioning, we found a 
rotation period P = 8.5 + 0. 1 days for the complete series as well 
as for both 2006 and 2007 seasons, with a confidence level larger 
than 99%. This is consistent with the period determination from 
ASAS data. The maximum light curve amplitude was 0.02 mag, 
whereas photometric accuracy was 0.01 mag. In Fig. 8 we show 
a summary of the period search on Super WASP data. We refer 
the reader to Messina et al. (2010) for a more detailed descrip- 
tion of the photometric data analysis methods. According to the 
criteria adopted in Messina et al. (2011) for assigning a confi- 
dence level to the rotation period, the period of HD 8049 can be 



http://www.wasp.le.ac.uk/public/ 



considered confirmed. The P - 8.3 days rotation period is com- 
patible with the projected rotational velocity (3.5 ±1.2 km/s, see 
Torres et al. 2006) for stellar radii R > 0.6 Rq. 

4.2. Age indicators 

Table 3 summarizes the measurements and corresponding ages 
as derived from the different indicators. From our analysis of 
the HARPS spectrum, we confirm that no lithium is detectable, 
as found by Torres et al. (2006). Two sets of log/?HK measure- 
ments are available: those from our CORALIE spectra and that 
by Gray et al. (2006). Within the uncertainties they give the 
same value (log7?HK = -4.25). The logarithm of the ratio of 
X-ray and bolometric luminosity logLx/i'boi was derived from 
ROSAT All-Sky Bright Source Catalog (Voges et al. 1999) as- 
suming the source IRXS JOl 1915.1-433803 is the X-ray coun- 
terpart of HD 8049 A. The X-ray emission should be dominated 
by the active K dwarf primary since HD 8049 B is cooler than 
WDs that show significant X-ray emission (Tg^ > 25000 K, 
Vennes 1999). The U, V and W heliocentric space velocities 
(given in a right handed coordinate system with U towards the 
galactic centre) are taken from Holmberg et al. (2009) adopting 
the trigonometric parallax of 29.79 ±1.21 mas (van Leeuwen 
2007). The uncertainty in the center of mass velocity due to or- 
bital motion spans from a minimum value of 200 m/s to a maxi- 
mum of 3 km/s, as determined from the simulations described in 
Sec. 5.1. 
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Fig. 7. Absolute magnitude in different bands (from top-left FUV, NUV, V and H) versus cooling age for white dwarf models of 
Vennes et al. (2011) (black lines) and sample a of nearby dwarfs collected by Giammichele et al. (2012). Stars represent objects 
with all magnitudes available and diamonds represent the others. Colors indicate the mass from the lowest (blue) to the highest 
(red ones). The continuous blue horizontal line represents the magnitude of HD 8049 B. The corresponding error bars are plotted as 
dashed blue lines. The plots show that the objects of the Giammichele et al. (2012) catalog are not peculiar and are well described 
by the theoretical models of Vennes et al. (201 1). 



The age estimates provided from these various indicators are 
quite discordant. Chromospheric HK emission, coronal emis- 
sion, and rotation period indicate an age between 100-400 Myr 
On the other hand, the lack of lithium indicates an age older than 
about 500 Myr. Additionally the kinematic parameters are quite 
far from the typical locus of young stars (Montes et al. 2001). 
The galactic orbit as given in Holmberg et al. (2009) has e = 0.19 
and Zma.K = 0.33 kpc. These properties are more typical of stars 
with age of about 3 Gyr, as estimated from the galactic orbits 
of stars with small age uncertainty in the Holmberg et al. (2009) 
catalog. These discrepancies strongly suggest that HD 8049 A is 
not a bona-fide young star but is rather a peculiar object. 

We identify the source of these anomalies to be the WD com- 
panion detected through direct imaging. The accretion of some 
amounts of mass and angular momentum by the central star at 
the time of mass loss from the WD progenitor might have spun- 
up the star, while increasing the level of magnetic activity. The 
presence of the WD cannot alter significantly the measurement 
of the rotation period and X-ray flux. Only minor effects on the 
logi^HK measurement are expected from the contamination of 
the spectrum. 

We conclude that the young rotation-activity age is not spu- 
rious. The central star should have been rejuvenated by the ac- 
cretion of some amount of mass and angular momentum at the 
time of mass loss from the WD progenitor This scenario is in- 
vestigated more thoroughly in Sec. 5.3. 



Table 3. Age indicators of HD 8049 A. The 2nd column Hsts the 
adopted value for the various parameters, and the 3rd column 
lists the corresponding reference; the 4th column contains the 
resulting age, and the last column has the adopted calibration for 
age determination. 



Indicator 


Measure 


Ref 


Age (Myr) 


Ref. 


Li EW (mA) 





1,2 


>500 


1 


log.RHK 


-4.25 ± 0.05 


1,3 


90 


4 


logLx/Lto, 


-4.24 


1 


182 


4 


Ao. (d) 


8.3 + 0.1 


1 


360 


4 


P,ot (d) 






380 ± 30 


5 


U, V, W (km/s) 


18,-47,-28 


6 


old (few Gyr) 


1 



References. 1 : this paper; 2: Torres et al. (2006); 3: Gray et al. (2006); 4: 
Mamajek & Hillenbrand (2008); 5: Delorme et al. (2011) 6: Holmberg 
et al. (2009) 



4.3. Mass determination 

Using our spectroscopic temperature and metallicity results, we 
employ the PARAM interface^ (da Silva et al. 2006) to conclude 



' http://stev.oapd.inaf.it/param 
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Fig. 6. Top: Observed magnitudes of the system plotted over the 
synthetic spectra. Blue points are those taken in consideration 
during the SED fitting of the WD. Triangle shaped ones are not 
considered during the PASTIS integration. Values of the magni- 
tudes of the integrated system are represented in red and the SED 
in black. The primary star SED is represented in green. Bottom: 
Two parameters joint posterior distributions for HD 8049 B. The 
68.3%, 95.5% and 99.7% confidence regions are denoted by 
light grey, dark gray and black areas respectively. Single PDF 
(Posterior Distribution Function) of each combination is shown 
along the sides. 



a mass of 0.804±0.023 M^. The mass of the convective envelope 
of the star, following MuiTay et al. (2001), is 0.086 Mq. 

^A. Abundance analysis 

It is well known that a considerable amount of mass can be 
transferred to a companion in close binary systems. This is the 
explanation for blue stragglers in clusters (McCrea 1964), Ba- 
stars and extrinsic C-stars (McClure et al. 1980; McClure 1984; 
McClure & Woodsworth 1990). Different signatures can be ob- 
served, depending on when and how much mass was transferred 
from the donor (in our case the white dwarf). The surface chem- 
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Fig. 8. Photometry time series of HD8049 from SuperWASP: 
left columns display time segments of magnitudes versus 
Heliocentric Julian Days (HJD) (the first panel only shows the 
complete series). Middle columns display the Lomb-Scargle pe- 
riodograms with indication of the rotation period. Right columns 
display the light curves phased with the rotation period and the 
first HJD as initial epoch. Solid curves represent the sinusoidal 
fit with the rotation period. The labels report the ratio of average 
residual from the fit and the rotational modulation amplitude. 



ical composition of the current primary can therefore be used to 
better constrain the history of this system. 

In order to better characterize the primary, and to explore 
possible signatures of accretion of material processed by the 
WD progenitor, we derived the stellar spectroscopic parameters 
and measured the abundances of several elements. To this aim, 
we used the HARPS high resolution spectrum of HD 8049 A 
(Sec. 2.2) and adopted the prescriptions of Biazzo et al. (2012). 
For a detailed description of the method, we refer to that paper 

The derivation of the iron abundance of our target was per- 
formed differentially with respect to the Sun. To this purpose, we 
analyzed a Ganymede spectrum acquired with HARPS at high 
S/N, obtaining logn(FeI)o = logn(FeII)0 = 7.53 ± 0.05. In the 
end, we found the following stellar parameters and iron abun- 
dance for HD 8049 A: effective temperature T^ff - 5050 ± 50 K, 
surface gravity log g - 4.50±0.10dex, microturbulence velocity 
^ =1.42 + 0.05 km/s, [Fel/H] = -0.10 + 0.02 dex, and [Fell/H] 
= -0.10 ± 0.05 dex. 

Looking for potential evidence of pollution by the ejecta of 
the precursor of the WD, we derived abundances of the s-process 
elements'^ Y, Ba, and La, which are mainly produced by low- 
mass AGB stars (~ 1.5-3 Mq, e.g., Busso et al. 1999). The abun- 
dances of Y, Ba, and La were derived as in D'Orazi et al. (2012). 
We also gathered Cu abundances, whose 70- 80% abundance 



** Elements heavier than iron are mostly synthesized via neutron- 
capture processes, that can be slow (s) and rapid (r), where slow and 
rapid is compared to the theyS-decay timescale (Burbidge et al. 1957). 
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in the Solar System is ascribed to massive stars (Kappeler et al. 
2011). We adopted the hyperfine splitting from Steffen (1985) 
and the isotopic mixture from Simmerer et al. (2003) (69% for 
^^Cu and 31% for^^Cu. We found [Y/Fe] = -0.15+0.10, [Ba/Fe] 
= -0.10+0.10, [La/Fe] = -0.12+0.10, and [Cu/Fe] = 0.00+0.07. 
Errors are those derived from the spectral synthesis. 

We also measured the abundances of Na and Al because 
these elements are expected to be overabundant in the ejecta 
of intermediate-mass AGB stars (Ventura et al. 2001), which is 
another potential precursor of HD 8049 B. We found [Na/Fe] - 
0.05 ± 0.04 and [Al/Fe] = 0.13 + 0.03, indicating no significant 
pollution from this kind of star. 

5. System properties 

5.1. Orbital properties 

To constrain the possible orbit of the object we performed a 
Monte Carlo simulation following the example of Desidera et al. 
(201 1). The simulation steps through 1.8x 10^ random orbital el- 
ements and rejects all the orbits that do not fit both the RV trend 
and the astrometric data. We adopted Ma - 0.804 ± 0.023 Mq 
and Mb = 0.56 ± 0.08 M^ as described, in Sec. 4.3 and 3.2. 
The simulation found ~28000 orbits consistent with the obser- 
vational data. We present the allowed orbital elements in Figure 
9. 

The distribution of each element is shown in the histograms 
of Figure 10. The distribution of the inclinations is bimodal and 
has its peaks on the values of 5 deg and 65 deg. The argument 
of periastron a> and the longitude of the node Q have the highest 
probability around values of 180 deg and 250 deg respectively. 
The most probable period is around 300 yr, corresponding to a 
semi-major axis of ~ 50 AU. Among the possible configurations 
of the orbit, there is a lack of edge-on inclinations that can pro- 
duce the correlation plots of w and Q.. This is consistent with 
the RV trend. The orbit is likely eccentric, with a most probable 
value e = 0.55. 

5.2. Broad band pliotometry 

Table 4 lists the available broad band photometry for the system. 
A very bright FUV magnitude is listed in the GALEX catalog. 
The offset between the optical and GALEX position of 1.67" 
along the SE direction indicates that the FUV source is better 
coincident with the WD rather than the K2 star. Such a bright 
FUV magnitude is not expected from the primary, even taking 
into account the star's high level of chromospheric activity (see 
e.g. Fig. 6 in Smith & Redenbaugh 2010). Corrections for chro- 
mospheric contribution of the primary have been calculated us- 
ing equations 2 and 4 of Findeisen et al. (2011); we assumed a 
B-V color for the K star as given in Table 15.7 of Allen (1977); 
we assumed log7?HK = -4.25. The final derived errors take into 
account the dispersion of the points used to derive the mentioned 
relations. 

Based on the composite U magnitude of the unresolved sys- 
tem, we determined an individual U magnitude following the 
procedures described in Sec. 2.3. H magnitude was obtained 
with NACO observations (see Sec. 2.1). 

5.3. System iiistory 

As mentioned in Sec. 4.2, the discrepant results of various age 
indicators for HD 8049 A might be qualitatively explained by the 
accretion of some amounts of mass and angular momentum from 



Table 4. Broad band photometry of HD 8049 



Band 


Mag 


Ref. 


A+B 


FUV 


12.645 ± 0.004 


1 


NUV 


13.454 + 0.003 


1 


FUV-NUV 


-0.809 ± 0.005 


1 


U 


10.098 + 0.031 


2 


B 


9.626 + 0.028 


2 


V 


8.734 + 0.027 


2 


B-V 


+0.892 + 0.007 


2 


V-I 


+0.95 ± 0.02 


3 


U-B 


+0.472 + 0.013 


2 


J 


7.077 ± 0.027 


4 


H 


6.649 + 0.059 


4 


K 


6.523 + 0.031 


4 


V-H 


+2.04 + 0.06 


2,4 


V-K 


+2.17 + 0.04 


2,4 


A 


FUV 


19.68 + 0.57 


5 


NUV 


15.12 + 0.22 


5 


U 


10.30 + 0.03 


5 


V 


8.747 + 0.030 


5 


H 


6.65 + 0.07 


5 


B 


FUV 


12.65 + 0.33 


5 


NUV 


13.72 + 0.21 


5 


u 


12.02 + 0.04 


5 


V 


13.5 + 0.5 


6 


H 


13.75 + 0.12 


5 



References. 1: GALEX, AB magnitude system (Martin et al. 2005) 2: 
Cousins (1983); 3: Ferryman et al. (1997); 4: 2MASS (Skrutskie et al. 
2006) ; 5: this paper; 6: van den Bos (1929); 



its WD companion HD 8049 B. Comparable cases have been dis- 
covered in the past. Jeffries & Stevens (1996) showed that the 
accretion of slow wind from the AGB progenitor of the WD can 
transfer sufficient angular momentum to speed up the companion 
even for orbital separation up to 100 AU. Cases of pairs formed 
by a WD and a fast rotating, barium-rich star have been reported 
by Jeffries & Smalley (1996). Another case is represented by 
the binary GJ86, which is known to host a giant planet and a 
WD companion (GJ86B; Els et al. 2001; Mugrauer & Neuhauser 
2005). GJ86 also shows anomalies in the age indicators that are 
similar, even if less extreme, than those we found for HD 8049; 
the GJ86 age indicators yield a younger age from X-ray and 
H&K than from kinematic (see Desidera & Barbieri 2007, App. 
B). Another more direct indication of wind accretion on a com- 
panion during the AGB phase is represented by the accretion 
disk around Mira B (Karovska et al. 1997; Ireland et al. 2007), 
which lies at a projected separation of about 70 AU, which is 
slightly larger than the most probable value for HD 8049. 

Following these works we can speculate that part of the ma- 
terial lost through stellar wind by the WD progenitor was ac- 
creted by HD 8049 A, increasing its rotation rate and mimick- 
ing a young age for indicators that are directly or indirectly tied 
to stellar rotation (rotation itself. X-ray and H&K emission). 
Specifically, we expect that the activity and rotation age of the 
system is related to the time elapsed since the accretion event, 
i.e. the cooling age of the WD. However, the activity-rotation 
age might be larger than the cooling age of the WD if the pri- 
mary was not spun-up to the fast rotational velocity typical of a 
very young star. In this scenario, increasing the mass of the WD 
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progenitor allows for larger amounts of accreted material, both 
because of large amounts of mass loss due to stellar wind at the 
end of AGB phase, and also because of the larger widening of the 
binary orbit due to mass loss. A more massive progenitor implies 
also a younger age for the system. On the other hand, the lack 
of significant abundance anomalies on HD 8049 A (Sect. 4.4) 
indicates that either the s-process enhancement of the material 
lost by the AGB progenitor is small, which may arise only for 
masses lower than 1.5 Mq, or that the amounts of accreted ma- 
terial is small with respect to the heavy element content of the 
convective zone. 

To explore the possible original configuration of the system, 
we used the BSE package (Hurley et al. 2002), evolving binary 
systems with different primary masses and orbital parameters. 
We used their adopted initial-final mass relationship for the WD 
and default parameters to take into account mass loss and wind 
accretion. These elements have significant uncertainty so the re- 
sults are intended only for a qualitative view of the possible con- 
figuration. 

A system with a 1.2-1.5 M© star initially at 30 AU will pro- 
duce a WD of roughly the correct mass with a semi-major axis 
close to the most probable value, for current orbits presented in 
Sect. 5.1. The resulting age of the system (3-6 Gyr) is compatible 
with the kinematic properties of HD 8049. About 0.005 Mq ai-e 
accreted by the other component in this case, which is about 5% 
of the mass of the stellar convective zone as derived following 
the prescriptions by Murray et al. (2001). 

We checked that the results using the more recent semi- 
empirical relationship by Salaris et al. (2009) are fully consis- 
tent with those derived using BSE. Salaris et al. (2009) found 
that WDs with masses in the range 0.56-0.64 Mq are found in 
NGC 7789 (estimated age of 1.5 Gyr) and older clusters while 
those in younger ones, like the Hyades and Praesepe (age of 
about 0.65 Gyr), are more massive. 



6. Conclusions 

The substellar candidate around HD 8049 A, initially suspected 
based on NACO //-band photometry and astrometry, turned out 
to be a WD companion rather that a BD. 

With direct imaging, second epoch observations of the ob- 
ject, we confirmed a bound companionship and determined the 
WD mass through its luminosity. 

The M/L relation for exoplanets strongly depends on the 
age of the system (see e.g. models from BarafFe et al. 2003; 
Madhusudhan et al. 201 1). For this reason, samples of stars scru- 
tinized for low mass companions through high contrast imag- 
ing usually include a high fraction of presumably young stars. 
However, a misinterpretation of the age of the host star may pro- 
duce false alarms and/or an erroneous calculation of the mass of 
the companion. This was the case of HD 8049 A, a star which 
appears rejuvenated by the accretion of significant angular mo- 
mentum from the ejecta of the progenitor of a WD. The faint 
companion was confused with a 35 Mjup brown dwarf when the 
age indicated by rotation and activity was considered. 

From radial velocity analysis we found a trend that was com- 
patible with a more massive object rather than a brown dwarf. 
Indeed, a comprehensive analysis of age indicators revealed in- 
consistencies. The final clue that unequivocally confirms the WD 
hypothesis is the observation of the system in U- band. At this 
wavelength, a bright object is observed at the same astrometric 
position of the NACO candidate. The ultraviolet magnitude of 
this object is 1 .7 mag fainter than the host star. 



We found parameters for the WD assuming an atmosphere 
made of pure hydrogen and fitting photometry observation 
points. Results are: T^ff = 18800 ± 2100 K, Mwd = 0.56 ± 
0.08 Mq, log^ = 7.86 ± 0.16 and age feooi = 63+^^ Myr 

We made a Monte Carlo simulation, following the example 
of Desidera et al. (201 1), to constrain the orbital parameters. We 
found that the WD moves in a likely eccentric orbit with a most 
probable period of 300 yr (~ 50 AU). 

In our scenario, the existence of the WD companion (with 
//-band magnitude compatible with a sub-stellar object) and the 
spurious young age of the star, as derived from rotational pe- 
riod, log/?HK and X-ray emission, are intrinsically linked. As 
discussed in this paper, a way to identify the spurious young stars 
that have been affected by wind accretion is to compare the ages 
resulting from several age indicators, paying attention to discrep- 
ancies between rotation-activity ages and lithium and kinematic 
ages. In the case that the sample is focused on young stars, such 
a false alarm is related to a WD with a short cooling age, which 
explains a young rotation-activity age for the companion. Such 
objects are most easily identified in the UV bandpass, as done in 
our study exploiting archive GALEX observations. 

Resolved systems with white dwarf companions close 
enough to have harboured some accretion phenomena are in- 
teresting targets for further constraining the wind accretion oc- 
curring in moderately wide binaries during the asymptotic giant 
branch (AGB) phase of the WD progenitor They also enable 
investigations of the maximum binary separation at which Ba- 
stars can be observed. Finally, our results confirm that there is a 
population of WD companions awaiting discovery even at small 
distance from the Sun, as suspected by Holberg (2009), and that 
AO imaging can play a relevant role for unveiling it. 

As nearby young stars are prime targets for direct imaging 
planet searches, and new instrumentation optimized to this sci- 
ence goal is about to start operations, one might ask what is the 
expected frequency of the kind of false alarm we have found 
in NACO-LR To this purpose, we may first consider the fre- 
quency of Ba-stars, which are systems similar to those consid- 
ered here, but with initially smaller separation and shorter pe- 
riod. Ba-stars are the result of mass transfer from an AGB com- 
panion (the original primary of the system) to the small mass 
star that is currently observed. For our considered scenarios, the 
system is observed at a late stage in its development, when the 
AGB star has evolved into a WD. From a survey of ~200 slightly 
evolved F dwarfs, Tomkin et al. (1989) discovered two barium 
dwarfs, and derived a rough frequency of ~1% Ba-stars among 
F dwarfs; Edvardsson et al. (1993) found 6 dwarf barium stars 
among 200 F stars (i.e. 3%), but their binary nature remains to 
be proven. North & Duquennoy (1991) found a frequency of 0.5 
to 1.0% Ba dwarfs among FV stars. McConnell & Gast (1972) 
estimated that 1 % of the G and K giants in the vicinity of the 
Sun are Ba-stars. From a 200 star sample, Williams (1975) de- 
rived a frequency of Ba-stars among G-K giants between 1 and 
4%. Bidelman (1985) obtained a frequency between 0.5 and 1%. 
Most Ba-stars are in binary systems with orbital periods typi- 
cally longer than 100 days but less than 10 years (McClure et al. 
1980; McClure 1984; McClure & Woodsworth 1990), though 
Jorissen et al. (1998) reported a few cases with periods outside 
this range. Orbits typically have low but non-zero eccentricity, 
linked to the fact that they are likely not produced by Roche 
lobe overflow; the origin is probably due to stellar wind. This is 
further supported by evolutionary models (Han et al. 1995). 

We may consider the period and secondary mass distribution 
of binary stars, and estimate the fraction of systems that are of 
potential interest. Integrating the distribution of binary separa- 
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tion by Duquennoy & Mayor (1991), it comes out that at least 
half of the binaries containing a WD have a separation large 
enough that the star does not become a Ba-star Hence, the frac- 
tion of G-K stars with a WD companion should be of the order of 
a few percent. This is consistent with arguments based on the ini- 
tial mass function: roughly a few percent of the G-K stars should 
have a companion originally more massive if we assume that the 
mass of the two stars is independent of each other, which is prob- 
ably a good approximation at least for solar type stars. Objects 
like HD 8049 A should therefore most likely be G-K stars with 
a WD companion and an original separation typically slightly 
larger than that required to become a Ba-star. They should also 
have experienced some mass and angular momentum transfer 
Hence there should not be large separations. 

Objects like HD 8049 A may very well pass through the pre- 
selection of the direct imaging candidates, failing to show up as 
spectroscopic binaries. Periods are long enough that RV varia- 
tions are very small. Furthermore, they will not show up in as- 
trometric binary catalogs, because periods are long that orbital 
motion is not distinguished from proper motion. Finally, they 
might also be overlooked as visual binaries because the con- 
trast is large. Inclusion in the sample of direct-imaging planet 
searches also depends on the adopted age limits. For surveys fo- 
cused on young stars, such as our NACO one, only companions 
of WDs with short cooling times might be selected. 
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Fig. 9. Monte Carlo simulation results showing orbital parameters for orbital solutions compatible with the observational data. 
Colors indicate different levels of ;t'-squared values (decreasing from black to blue and finally red). The plots show that some 
parameters are well constrained in a region of parameters space while others can span all possible values. See 5.1 for information 
on the Monte Carlo procedures. 
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Fig. 10. Histograms of orbital parameters derived with a Monte Carlo simulation following the example of Desidera et al. (201 1). 
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